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ABSTRACT 

We investigate five different models for the dark matter halo bias, ie., the ratio of 
the fluctuations of mass tracers to those of the underlying mass, by comparing their 
cosmological evolution using optical QSO and galaxy bias data at different redshifts, 
consistently scaled to the WMAP7 cosmology. Under the assumption that each halo 
hosts one extragalactic mass tracer, we use a x 2 minimization procedure to determine 
the free parameters of the bias models as well as to statistically quantify their ability 
to represent the observational data. Using the Akaike information criterion we find 
that the model that represents best the observational data is the Basilakos & Plionis 
(2001; 2003) model with the tracer merger extension of Basilakos, Plionis & Ragone- 
Figucroa (2008) model. The only other statistically equivalent model, as indicated 
by the same criterion, is the Tinker et al. (2010) model. Finally, we find an average, 
over the different models, dark matter halo mass that hosts optical QSOs of: Mh — 
2.7(±0.6) x 1O 12 /i _1 M0, while the corresponding value for optical galaxies is: Mh — 
6.3(±2.1) x IO^^Mq. 



INTRODUCTION 



It is of paramount importance for cosmological and galaxy 
formation studies the understanding of how galaxies and 
other extragalactic mass-tracers relate to the underlying dis- 
tribution of matter. The current galaxy formation paradigm 
assumes that galaxies form within dark matter haloes, iden- 
tified as high-peaks of an underlying initially Gaussian den- 
sity fluctuation field, and that they trace in a biased man- 
ner such a field (eg., Kaiser 1984; Bardeen et al. 1986). A 
formation process of this sort can explain the difference of 
the clustering amplitude between the different extragalactic 
mass tracers (galaxies, groups and clusters of galaxies, AGN, 
etc) as being due to the different bias among the underlying 
density field and that of the dark matter (DM) halos that 
host the mass tracers. 

In order to quantify such a difference, one can use the so- 
called linear bias parameter b, which for continuous density 
fields is defined as the ratio of the fluctuations of the mass 
tracer (5 tr ) to those of the underlying mass (<5 m ): 

b=^, (1) 
Based on this definition one can write the bias parameter in 



a number of equivalent ways: (a) as the square root of the 
ratio of the two-point correlation function of the tracers to 
the underlying mass: 

since £(r) = (<5(x)<5(x + r)}, in which case one considers 
the large-scale correlation function (ie., scales ^ lh^ 1 Mpc), 
corresponding roughly to the so-called halo-halo term of the 
DM halo correlation function (eg., Hamana, Yoshida, Suto 
2002), and (b) as the ratio of the variances of the tracer 
and underlying mass density fields, smoothed at some linear 
scale traditionally taken to be 8 Mpc (at which scale 
the variance is of order unity): 



b = 



C"8,tr 
08, m 



CD 



since af = £(0) = (<5 2 (x)). 

A further important ingredient in theories of structure 
formation, is the cosmological evolution of the DM halo bias 
parameter (eg., Mo & White 1996; Tegmark & Peebles 1998, 
etc). A large number of such bias evolution models have 
been presented in the literature and the aim of this work is 
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to compare them using as a criterion how well do they fit 
the observed bias, at different redshifts, of optical QSOs and 
galaxies. In such a comparison we will make the simplified 
assumption that each DM halo hosts one mass tracer. This 
is consistent with the definition of the linear bias, where one 
uses either the large-scale correlation function (which cor- 
responds to the halo-halo term) or the smoothed to linear 
scales variance of the fluctuation field, while any residual 
non-linearities will probably be suppressed in the ratio of 
the tracer to underlying mass correlation functions or vari- 
ances. Further suppression of non-linearities, introduced for 
example by redshift-space distortions, can be achived using 
the integrated correlation function within some spatial scale; 
see discussion in section 2 below. 

There are two basic families of analytic bias evo- 
lution models. The first, called the galaxy merging bias 
model, utilizes the halo mass function and is based on 
the Press-Schechter (1974) formalism, the peak-background 
split (Bardeen et al. 1986) and the spherical collapse model 
(Cole & Kaiser 1989; Mo & White 1996, Matarrese et 
al. 1997; Moscardini et al. 1998; Sheth & Tormen 1999; 
Valageas 2009; 2011). Cole & Kaiser (1989) found for the 
bias evolution that 

b(M, z) = 1/(1 + z) - 1/ [1.68(1 + z)] + 1.68(1 + z)/o 2 {M) , 

where a 2 (M) is the variance of the mass fluctuation field, 
while Mo & White (1996) derived for an Einstein-de Sitter 
universe that 

b(z) = 0.41 + [6(0) - 0.41] (1 + zf . 

Mo, Jing & White (1997) extended the previous study in 
the quasi-linear regime by taking into account high order 
correlations of peaks and halos. Similarly, Matarrese et. al. 
(1997) estimated the bias in a merging model where the 
halo mass exceeds a certain threshold. They found that for 
an Einstein-de Sitter universe: 

b(z) = 0.41 + [6(0) - 0.41] (1 + zf , 

while Moscardini et. al. (1998) generalized the above bias 
evolution model for a variety cosmological models. 

Many studies have compared the prediction of the merg- 
ing bias model with numerical simulations and beyond an 
overall good agreement, differences have been found in the 
details of the halo bias. For example, the spherical collapse 
model under-predicts the halo bias for low mass halos and 
fails to reproduce the dark matter halo mass function found 
in simulations. To solve this problem, Sheth, Mo & Tor- 
men (2001) extended their original model to include the 
effects of ellipsoidal collapse. However according to Tinker 
et. al. (2010), this model under-predicts the clustering of 
high-peaks halos while over-predicts the bias of low mass 
objects. Furthermore, Manera et al. (2009) and Manera & 
Gaztanaga (2011) find that the clustering of massive halos 
cannot be reproduced from their bias calculated using the 
peak-background split. 

Such and other differences have lead to other modi- 
fications of the models, either suggesting new fitting bias 
model parameters (eg., Jing 1998; Tinker et al. 2005), or 
new forms of the bias model fitting function (eg. Seljak & 
Warren 2004; Pillepich et al. 2010; Tinker et al. 2010) or 
even a non-Markovian extension of the excursion set theory 
(Ma et al. 2011). A further step was provided by de Simone, 



Maggiore & Riotto (2011), who incorporated the effects of 
ellipsoidal collapse to the original Ma et al. model, which is 
based on spherical collapse. 

The second family of bias evolution models assumes 
a continuous mass-tracer fluctuation field, proportional to 
that of the underlying mass, and the tracers act as "test 
particles". In this context, the hydrodynamic equations of 
motion and linear perturbation theory are applied. This fam- 
ily of models can be divided into two sub-families: 
(a) The so-called galaxy conserving bias model uses the con- 
tinuity equation and the assumption that tracers and un- 
derlying mass share the same velocity field (Nusser & Davis 
1994; Fry 1996; Tegmark & Peebles 1998; Hui & Parfey 2007; 
Schaefer, Douspis & Aghanim 2009). Then the bias evolution 
is given as the solution of a 1st order differential equation, 
and Tegmark & Peebles (1998) derived: 



b(z) = 1 + [6(0) - 1]/D(z) , 



where 6(0) is the bias factor at the present time and D(z) the 
growing mode of density perturbations. However, this bias 
model suffers from two fundamental problems: the unbiased 
problem ie., the fact that an unbiased set of tracers at the 
current epoch remains always unbiased in the past, and the 
low redshift problem ie., the fact that this model represents 
correctly the bias evolution only at relatively low redshifts 
z< 0.5 (Bagla 1998). Note that Simon (2005) has extended 
this model to also include an evolving mass tracer population 
in a ACDM cosmology. 

(b) An extension of the previous model, based on the ba- 
sic differential equation for the evolution of linear density 
perturbations, which implicitly uses that mass tracers and 
underlying mass share the same gravity field, and on the 
assumptions of linear and scale-independent bias, provides 
a second order differential equation for the bias. Its approx- 
imate solution provides the functional form for the cosmo- 
logical evolution of bias (Basilakos & Plionis 2001; 2003 and 
Basilakos, Plionis & Ragone-Figueroa 2008; hereafter BPR 
model) . The provided solution applies to cosmological mod- 
els, within the framework of general relativity, with a dark 
energy equation of state parameter being independent of cos- 
mic time (ie., quintessence or phantom). An extension of this 
model to engulf also time-dependent dark energy equation of 
state models, including modified gravity models (geometric 
dark energy), was recently presented in Basilakos, Plionis & 
Pouri (2011). 

The outline of this paper is as follows. In section 2 we 
present the data that we will use, we review the basic tech- 
niques used in measuring the bias from samples of extra- 
galactic objects and we will present the rescaling method 
used in order to transform different bias data to the same 
(WMAP7) cosmology (ie., flat ACDM with Q m = 0.273 and 
g 8 = 0.81). In section 3 we introduce the different bias evo- 
lution theoretical models that we will investigate, while in 
section 4 we present our results and discussion. The main 
conclusions are presented in section 5. In the Appendix we 
discuss the simulations used to fit the free parameters of the 
BPR model, as well as the cosmological dependence of these 
parameters. 
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2 MASS TRACER BIAS DATA 

The mass tracers that we will use in this work are op- 
tical QSOs and galaxies, for which their linear bias with 
respect to the underlying mass is available as a function of 
redshift. In particular, we will use: 

• The 2dF-based QSO results of Croom et al. (2005), 
which are based on spectroscopic data of over 20000 QSOs 
covering the redshift range 0.3 ^ z ^ 2.2 and on a ACDM 
cosmology with Q m = 0.27 and erg = 0.84. 

• The SDSS (DR5) QSO (z< 2.2) results of Ross et 
al.(2009) based on spectroscopic data of ~30000 QSOs and 
and on a ACDM cosmology with fi m = 0.237 and <jg = 
0.756. 

• the SDSS (DR5) QSO results of Shen et al. (2009), who 
used a homogeneous sample of ~38000 QSOs within 0.1 ^ 
2^5 and on a ACDM cosmology with fi m = 0.26 and 
erg = 0.78. In this case we will use only their z~^ 2.2 results 
to avoid including in our analysis correlated measurements 
of the bias, for the redshift range covered also by the Ross 
et al. analysis. 

Although there are other QSO bias data available, like the 
Myers et. al. (2006) analysis of 300000 photometrically clas- 
sified SDSS DR4 QSOs, within 0.75 ^ z sC 2.8, we do not 
include them in our analysis in order to avoid, in the redshift 
range studied, as much as possible correlated bias measure- 
ments. 

As far as galaxy data are concerned, we will use the 
bias results of Marinoni et. al. (2005), which are based on 
3448 galaxies from the VIMOS-VLT Deep Survey (VVDS), 
cover the redshift range: 0.4 ^ z ^ 1.5 and use a ACDM 
cosmology with Q. m = 0.3 and erg = 0.9. 

Although in the next section we sketch the usual pro- 
cedures used to estimate the linear bias of a sample of ex- 
tragalactic mass tracers, we would like to stress that for the 
QSO data used in this work, the corresponding authors, in 
order to minimize non-linear effects, have estimated the in- 
tegrated correlation function for scales > l/i -1 Mpc, which 
in the usual jargon corresponds roughly to the halo-halo 
term of the DM halo correlation function. As for the VVTS 
galaxy bias data, Marinoni et al., devised a procedure to 
estimate the bias of a smooth galaxy density field in pencil 
beam surveys, disentangling the non-linear effects, and thus 
the bias values used in this work correspond to the linear 
bias. 



2.1 Estimating the tracer bias at different 
redshifts 

Although we will use the bias data provided by the previous 
references, for completion we briefly present here the basic 
methodology used to estimate the bias of some extragalactic 
mass tracer at a redshift interval z ± Sz, using any of the 
basic definitions of eq.(l)-(3). 

The first issue that one has to keep in mind is that 
what we measure from redshift catalogues is the redshift- 
space distorted value of either the tracer correlation func- 
tion, £tr(s), or the variance of the tracer density field crf it r,s 
(the index s indicates redshift-space distorted spatial sepa- 
rations, while the index r indicates true spatial separations). 
One needs to correct for such distortions, resulting from the 



peculiar velocities of the mass tracers, in order to recover the 
true spatial value of either measures. Kaiser (1987) provides 
such a correction procedure which entails in dividing the di- 
rectly measured from the data tracer correlation function or 
variance with a function F({l m , f^A, b, z), given by (see also 
Hamilton 1998 and Marinoni et al. 2005): 

F(n m ,n A ,b,z) = i + |/?( 2 ) + I/3 2 ( z ) (4) 

with p(z) = W m {z)/b(z), and 7 = 6/11 for the ACDM (eg., 
Wang & Steinhardt 1998; Linder 2005), which implies that 
P(z) = n^Eiz)- 12 ' 11 ^ + z^'^/biz). Therefore the re- 
lation between the redshift-space and real-space measures 
used to estimate the bias parameter is: 

c^ = ^>) = ( m ' A ' ' } ( ) 

Then combining equations (J2J) or (J3| with (T5J) and $5$ 

provides the real-space bias factor according to: 



b(z) = 



ftr(s,«) 


^ 12/11 {z) 


1/2 


n™ /n (*) 


_£ m (r, z) 


45 




3 




4^ 2/11 (*)~ 


1/2 


nT\z) 




45 




3 



where fm(r) and erf m are the corresponding correlation 
function and variance of the underlying dark matter distri- 
bution, given by the Fourier transform of the spatial power 
spectrum P(k) of the matter fluctuations, linearly extrapo- 
lated to the present epoch: 

and 

7-j2 / \ poo 

<rl m (z) = J k 2 P(k)W 2 (kR s )dk , (8) 

with D(z) the normalized perturbation's growing mode (ie., 
such that -D(O) = 1), P(k) the CDM power spectrum given 
by: 

P(k) = P k n T 2 (tl m ,k) , (9) 

with T(Q m ,k) being the CDM transfer function (Bardeen 
et al. 1986; Sugiyama 1995; Eisentein & Hu 1998), n the 
slope of the primordial power-spectrum (which according to 
WMAP7 is = 0.967) and W(kR 8 ) the Fourier transform of 
the top-hat smoothing kernel of radius R = Rg = 8/t -1 Mpc, 
given by W(kR 8 ) = 3(sinfcJ?g - fci? 8 cosfcflg)/(fc7? s ) 3 . 

Now, although in the case of using eq.(3), the erg vari- 
ance is free of non-linear effects by definition, this is not so 
when using the correlation function approach (eq. 2). There- 
fore, in order to minimize nonlinear effects at small separa- 
tions one can replace £tr(s) in eq.© with the integrated 
correlation function, £tr(s). 

An alternative approach in order to avoid redshift-space 
distortions is to resolve the redshift-space separation, s, into 
two components, one perpendicular (r p ) and one parallel 
(n) to the line-of-sight (see Davis & Peebles 1983) and then 
estimating the 2-point projected correlation function w p (r p ) 
along the perpendicular dimension (within some range of the 
parallel dimension, say TT m i n < n < Tv ma x), which is related 
to the spatial correlation function, £(r), according to: 
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w p (r p ) 



£(r p ,7r)d7r 



T rair , \T 



-dr 



(10) 



where ir — \Ad\ and r p — AdtanS/2, with Ad the radial 
comoving distance separation of any pair of mass tracers 
and 6 is angular separation on the sky of the pair members. 
As before, one can use the integrated correlation function, 
£, in order to minimize nonlinear effects. 

Additionally, one can also use the angular two-point 
correlation function, w(0), instead of £(s) or w p (r p ), in order 
to obtain £(r) via Limber's inversion, a procedure which 
also avoids the peculiar velocity distortions, but is hampered 
by the necessity of a priori knowing the redshift selection 
function of the mass tracers. 



2.2 Scaling the bias data to the same Cosmology 

Since different authors have estimated the optical QSO 
and galaxy bias using different cosmologies, we need to con- 
vert them to the same cosmological background in order to 
be able to use them consistently. As such we choose the re- 
cent WMAP7 cosmology (Komatsu et al. 2011). 

The procedure that we will follow uses the different erg 
power-spectrum normalizations (eq. 3). We wish to translate 
the value of bias from one cosmological model, say B, to 
another, say A. The definition of bias at a redshift z for 
these two different cosmologies are given by: 



b A (z) = 



and 



b B (z) 



CS,m,s(z) 



(11) 



(12) 



where the numerator is the real space value of as(z) esti- 
mated directly from the data, using also eq.© to correct 
for redshift space distortions. Dividing now equation (|11[) 
by ()12p , taking into account eq.([5]), and making the fair as- 
sumption that: 



C8,tr, s,a(z) ~ CTs.tr, s,b(^) 



(13) 



since the different cosmologies enter only weakly in the ob- 
servational determination of ffs.tr, through the definition of 
distances, we then have: 



b A (z) ~ b B (z) 



CS.m,B(z) 



0"8,n 



m.B, ^A,B 



,b B ,z) 



-F(fim,A,fiA,A,&A, Z) 



1/2 



(14) 



As it can be realized the required rescaled real-space bias, 
b A , enters also in the right hand side of the above equa- 
tion, making it rather complicated to analytically derive the 
full expression (using for example eq(S}. However, noting 
that the redshift-space distortion correction enters in the 
scaling of the bias, from one cosmology to another, as the 
square-root of the ratio of the F functions, the expected de- 
viation by using in the right-hand side of eq. (|14[) the crude 
approximation &a — bs, does not affect significantly this 
correction. In any case, the magnitude of the relevant cor- 
rection, (Fb I Fa) 1 ^ 2 , is extremely small, typically: ~ 0.8% 
at z — 0.24 dropping to ~ 0.1% at z — 2.1, and the overall 
scaling of the bias to different cosmologies is dominated by 
the ratio of the corresponding as(z) variances. 

We can facilitate our scaling procedure by using the 



as (z — 0) power-spectrum normalizations of the different 
models, a value always provided by the different authors. 
We therefore translate the values of os(z) to that at z — 
by using the linear growing mode of perturbations according 
to: as(z) = as(0)D(z). The final scaling relation from the B 
cosmology to that of A, therefore becomes: 



&a(z) — 6s(«) 



o-8,m,a(0) Db(z) 

O"8,m,A(0) Da(z) 



Ma 



,b B ,z) 



F(Q m , A , QA,A,b B , z) 



.(15) 



3 THEORETICAL BIAS MODELS 

Here we briefly describe the bias evolution models that we 
are going to compare. As discussed in the introduction, we 
separate the models in two families. The galaxy merging 
model family, based on the Press-Schether formalism and 
the peak-background split. The models that we will inves- 
tigate, representing this family, are the Sheth, Mo & Tor- 
men (2001) extension of the original Sheth & Tormen (1998) 
model (hereafter SMT), the Jing (1998) model, the Tinker 
et. al. (2010) (hereafter TRK) and the Ma et. al. (2011) 
model (hereafter MMRZ). All these models provide the bias 
of halos as a function of the peak-height parameter, v, where 

u = 6 c (z)/a(M h ,z) , (16) 

with Mh the halo mass, cr 2 (M/ l , z) the variance of the mass 
fluctuation field at redshift z, and S c (z) the critical linear 
overdensity for spherical collapse, which has a weak redshift 
dependence (see eq.18 of Weinberg & Kamionkowski 2003). 

The basic free parameter of these bias models, to be 
fitted by the data (although depending on the model one 
more parameter may be allowed to vary - see below), is v 
and through the evolution of o~(Mh,z) we will be able to 
derive the predicted bias redshift evolution, as well as the 
value of Mh ■ The latter value will be estimated by using the 
definition of a, eqs.Q and from which we have that: 



a(M h 



0-8 



r o °° dkk n+2 T 2 (k)W 2 (kR) 
/J 50 dkk"+ 2 T 2 (k)W 2 (kR a ) ' 



(17) 



with R = (3M h /4:TTp) 1/3 , R 8 = 8 h' 1 Mpc and p = 2.78 x 
10 11 O m ^ 2 M Q /Mpc 3 . 

The second family contains the so-called galaxy conserv- 
ing models and their extensions. These models are based on 
the hydrodynamical equations of motion and linear pertur- 
bation theory while the most general such model, that we 
will investigate, is that of Basilakos & Plionis (2001; 2003), 
extended to included a correction for halo merging in Basi- 
lakos, Plionis & Ragone-Figuera (2008). 

Below we present the functional form of the bias evolu- 
tion for each of the models that we will investigate: 

3.1 BPR 

Basilakos and Plionis (2001; 2003) using linear pertur- 
bation theory and the Friedmann-Lemaitre solutions derived 
a second-order differential equation for the evolution of bias, 
assuming that the mass-tracer population is conserved in 
time and that the tracer and the underlying mass share the 
same dynamics. 

The solution of their differential equation, for a flat cos- 
mology, was found to be (Basilakos & Plionis 2001): 
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b(z) = C 1 E(z) + C 2 E(z)I(z) + l (18) 
where E(z) = [fi m (l + z) 3 + Qa] 1/2 and 



I(z) = 



(19) 



The constants of integration depend on the halo mass, as 
shown in BPR, and they are given by: 



Ci(M h ) wai 



C* 2 (M h ) 




(20) 



(21) 



and the values of ai,2 and /3i,2 where estimated originally 
from a ~WMAP1 ACDM numerical simulation in BPR. We 
have since run a WMAP7 ACDM simulation, the details 
of which can be found in Appendix Al, and from which 
we have determined the new values of the ai t 2 and 
parameters (see Table Al). The cosmological dependence of 
these parameters is also discussed in Appendix A2. 

In BPR it was found that the original Basilakos & Plio- 
nis model could well reproduce the bias evolution for z < 3, 
but not at higher redshifts, indicating the necessity to extend 
the model to include the contribution of an evolving mass- 
tracer population. Such an extension was presented in BPR 
and it was based on a phenomenological approach, although 
the functional form for the effects of merging was based on 
physically motivated arguments (see Appendix A2 of BPR). 
To this end they introduced to the continuity equation an 
additional time-dependent term, ^(t), associated with the 
effects of merging of the mass tracers, which depends on the 
tracer number density, its logarithmic derivative and on 5tr- 
They parameterized this term using a standard evolutionary 
form: 



V(z) = AH {l + zY 



(22) 



where fj, and A are positive parameters which engulf the 
(unknown) physics of galaxy merging. The bias evolution is 
now given by: 

b(z) = CiE(z) + C 2 E(z)I(z) + y p (z) + 1 (23) 
where the additional halo-merging factor, y P {z), is given by: 

y p (z)=E(z) r(x)I(x)dx- E(z)I(z) r(x)dx (24) 

JO JO 

with t(z) = f(z)E 2 (z)/(l + zf and f(z) = A(fi - 2)(1 + 
z) ll E(z)/D(z). The values of both A and \x have been fitted 
using ACDM numerical simulations (see BPR) and it was 
found that \i ~ 2.5 — 2.6 independent of the halo mass, while 
A increases with decreasing halo mass, with A ~ 0.006 and 
for intermediate (ie., 10 13 < M h < 10 13 8 /i _1 M©) and higher 
mass halos, respectively. Evidently, the bias factor at z — 
is provided by: 



b(z) = Ci + C 2 I(z) + 1 



(25) 



Therefore in the current analysis we will leave Mh as a 
free parameter to be fitted by the data (BPR model) but we 
will also allow (a) the parameter ai to be fitted by the data, 
keeping A equal to its simulation based value (A = 0.006, 
BPR-I model), as well as the parameter A to be fitted by the 



data keeping ot\ equal to its simulation based value (ai = 
4.53, BPR-II model). 



3.2 SMT 

In Sheth et. al. (2001) the original work of Sheth & Tor- 
men (1999) was extended for the case of an ellipsoidal, rather 
than a spherical collapse. This new ingredient reduces the 
difference between theoretical expectations and simulation 
DM halo data. Considering ellipsoidal collapse the density 
threshold required for collapse, contrary to the spherical col- 
lapse case, depends on the mass of the final object. 

Using the ratio of the halo power spectrum to that of 
the underlying mass, they derived the functional form for 
the bias as: 



b{y) = 1 + 



yfa5 c (z) 
with f(u) = 



\yfa{au 2 ) + y/abiau 2 ) 1 c — f(v)] 
(au 2 Y 



{av 2 Y + 6(1 - c)(l - c/2) 



(26) 



where the free parameters where evaluated using N-body 
simulations to have values: a — 0.707,6 = 0.5 and c = 0.6. 
In particular the value of a was found to depend mostly on 
how the simulation DM halos were identified. In the case of 
a Friends of Friends (FoF) algorithm the value a — 0.707 
corresponds to the standard linking length of 0.2 times the 
mean inter-particle separation. Decreasing the linking length 
would increase the value of a and vice-versa (see discussion 
in SMT). Therefore, beyond the value of the DM halo mass, 
Mh (which will be estimated from the resulting value of 
u(Mh) via eq. (|17[) . we will also allow the parameter a to be 
fitted by the data. 



3.3 JING 

Jing (1998) used the clustering of simulation DM halos to 
derive an expression for the bias which is independent of the 
shape of the initial power-spectrum, being CDM or power- 
law. His corresponding expression is: 



b(u) 



05 



(0.06-0.02n) 



1 + 



(27) 



where n is the linear power spectrum index at the halo scale 
(ie., n = dhiP(k)/dlnk ~ -2 for M h ~ 10 13 /i _1 M Q ). The 
only free parameter of this model, to be fitted by the data, is 
the halo mass, Mh (which will be estimated from the fitted 
value of a(Mh) via eq |17|) . 



3.4 TRK 

Tinker et. al. (2010) measure the clustering of dark matter 
halos based on a large series of collisionless N-body simu- 
lations of the ACDM cosmology. DM halos were identified 
using the spherical overdensity algorithm by which halos are 
considered as isolated peaks in the density field such that 
the mean density is A times the density of the background. 
Their bias fitting function reads as: 



b{v) = 1 - A 



+ Bv" + Cv c 



(28) 
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where y = log 10 A. For the WMAP7 ACDM model the value 
which corresponds to the virialization limit is Aacdm — 
355. The rest of the parameters of the model are: A = 
1 + 0.24yexp[-(4/y) 4 ], B = 0.183, C = 0.019 + 0.107?/ + 
0.19exp[-(4/y) 4 ], a = 0.44y - 0.88, b = 1.5, c = 2.4. 

Therefore, we will fit the observational data using as a 
single free parameter the DM halo mass (Mh, derived via 
cr(Mh) in eq fTT)) and using with y — log 10 (AACDM). How- 
ever, we will also allow the latter parameter to be fitted by 
the data, simultaneously with Mh- 

3.5 MMRZ 

Ma et. al. (2011) extended the original Press-Schether ap- 
proach incorporating a non-Markovian extension with a 
stochastic barrier, where they assume that the critical value 
for spherical collapse is itself a stochastic variable, whose 
scatter reflects a number of complicated aspects of the 
underlying dynamics. Their model contains two parame- 
ters: k, which parameterizes the degree of non-Markovianity 
and whose exact value depends on the shape of the filter 
function used to smooth the density field, and a, the so- 
called diffusion coefficient, which parameterizes the degree 
of stochasticity of the barrier. Taking into account the non- 
Markovianity and the stochasticity of the barrier, the bias 
takes the form: 



b(u) 



1 + 



av — 1 + — 



yv 2 /2 



r(o,s§-; 



^6 C [1 - an + e Q " 2 /2r(0, ^)] 



(29) 



where a = (1 + Db) 1 , with Db the diffusion coefficient, 
and F(0, x) the incomplete gamma function. Without the 
stochasticity of the barrier one has Db = — ► a = 1. 

Ma et al. (2011) have found using N-body simulations 
that using a = 0.818 and k = 0.23 they can reproduce to 
a good extent both the simulation bias and the halo mass- 
function as a function of v. We will therefore use these pa- 
rameter values to fit the observational bias data in order to 
constrain Mh- Additionally, we will allow both a and Mh to 
be fitted simultanesouly by the data, using n — 0.44, since 
this is the value for a top-hat smoothing kernel in coordinate 
space. Note that the value of k appears to be almost inde- 
pendent of cosmology, as discussed in Maggiore & Riotto 
(2010). 



4 FITTING MODELS TO THE DATA 

In order quantify the free parameters of the DM halo 
bias models we perform a standard % 2 minimization pro- 
cedure between N bias data measurements, bi(z), with the 
bias values predicted by the models at the corresponding 
redshifts, b(p,z). The vector p represents the free parame- 
ters of the model and depending on the model their number 
is one or two. This procedure makes the simplistic assump- 
tion that each DM halo hosts one mass tracer, an assump- 
tion which is justified from the way the QSO and galaxy 
bias data have been estimated (see discussion in section 2). 

The x 2 function is defined as: 



E 



bi(z) - b(p,z) 



with (z) is the observed bias uncertainty. We have in total 
N = 22 measured bias data for the optical QSOs, spanning 
from z = 0.24 to z = 4, and N — 5 for the optical galaxies, 
spanning from z — 0.55 to z = 1.4. 

Note that the uncertainty of the fitted parameters will 
be estimated, in the case of more than one such parame- 
ter, by marginalizing one with respect to the other. How- 
ever, since such a procedure may hide possible degeneracies 
between parameters, we will also present the 1, 2 and 3<r 
likelihood contours in the parameter plane. 

Furthermore, since we will attempt to compare the dif- 
ferent models among them, the x 2 test alone is not sufficient 
for such a task, since different models may have a different 
number of free parameters. Instead we will use information 
criteria to compare the strengths of the different models, ac- 
cording to the work of Liddle (2004), a procedure that favors 
those models that give a similarly good fit to the data but 
with fewer free parameters (see for example Saini et al. 2004; 
Godlowski & Szydlowski 2005; Davis et al. 2007 and refer- 
ences therein). To this end we will use, the relevant to our 
case, corrected Akaike information criterion for small sample 
size (AIC C ; Akaike 1974, Sugiura 1978), defined, for the case 
of Gaussian errors, as: 



AIC C 



X + 2k + 2k(k -l)/(N-k- 1) 



(31) 



(30) 



where k is the number of free parameters, and thus when 
k = 1 then AIC C = Xmin + 2. A smaller value of AIC C indi- 
cates a better model-data fit. However, small differences in 
AIC C are not necessarily significant and therefore, in order to 
assess, the effectiveness of the different models in reproduc- 
ing the data, one has to investigate the model pair difference 
AAIC C = AICca-AICc,*,. The higher the value of |AAIC C |, 
the higher the evidence against the model with higher value 
of AIC C , with a difference |AAIC C |~ 2 indicating a positive 
such evidence and |AAIC C |<) 6 indicating a strong such ev- 
idence, while a value ^ 2 indicates consistency among the 
two comparison models. 



4.1 Optical QSO Results 

Here we fit the different bias evolution models to the scaled 
to the WMAP7 cosmology optical QSO data, described in 
section 2. It is important to note that all the bias models 
used in this work (except the BPR) have been studied as a 
function of the threshold v, eq. (|16[) . ie., in effect as a function 
of the variance of the fluctuation field and thus as a func- 
tion of halo mass, while the free parameters of most models 
have been fitted using z = simulations. In these models 
the redshift dependence of the bias comes mostly from the 
redshift dependence of the peak-height, v (see eq |16ll . 

We will present separately the results of the models with 
one free parameter, the halo mass, and the models with an 
additional free parameter, as discussed in the theoretical 
model presentation sections. 



4-1.1 One free parameter models 

In Table 1 we present the best fit model parameters based 
on the x 2 minimization procedure, with the first and sec- 
ond columns listing the fitted halo mass, Mh, derived using 
eq. (|17|l and the value of bias at z = 0, respectively. We 
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Model 




6(0) 


X 2 - /df 


AIC C 


BPR 


3.0 ±0.4 


1.02 


12.88/21 


14.88 


SMT 


3.2 ± 0.4 


1.07 


23.21/21 


25.21 


JING 


2.1 ± 0.2 


0.98 


18.00/21 


20.00 


TRK 


3.0 ±0.4 


1.00 


15.79/21 


17.79 


MMRZ 


2.2 ± 0.2 


0.87 


20.14/21 


22.14 



Table 1. Results of the x 2 minimization procedure between the 
optical QSO data (TV = 22) and bias models with one free pa- 
rameter (k = 1). 
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Figure 1. Main Panel: Comparison of the QSO bias data (open 
circles correspond to Croom et. al. 2005; filled circles to Ross et 
al. 2009 and filled triangles to the high redshift data of Shen et. al. 
2010), with the one free parameter bias model fits (BPR: black 
continuous line; SMT: magenta dashed line; JING: green long- 
short dashed line; TRK: red short-dashed line; MMRZ: blue dot- 
dashed line). Lower Panel: The relative difference between the 
BPR model and all the rest, A&bpr( 2 )- Inset Panel: The value 
°f X 2 — X 2 n i n as a function of halo mass, M/,, for the indicated 
bias models. 

also present the goodness of fit statistics, as discussed pre- 
viously (reduced x 2 and AIC C ). In the main panel of Figure 
1 we present the bias evolution models (different lines), us- 
ing the best fit parameter values listed in Table 1 together 
with the WMAP7-scaled optical QSO bias data. The inset 
panel of Fig.(l) shows that the resulting Mh values cluster 
around two, relatively similar, values: ~ 3 x 10 12 hT 1 M© and 
~ 2 x lO 12 /i _1 M0. In the lower panel we present the rela- 
tive difference between the BPR model and all the rest, ie., 
A6 B pr(2) = [bi{z) - &bpr(.z)]/6 B pr(,z). 

Some basic conclusions that become evident, inspecting 
also Table 1, are: 

• Although all one free parameter bias models appear 
to fit at a statistically acceptable level the optical QSO 
bias data, by far the best model is the BPR, which is the 
only model fitting also the highest redshifts (z£ 3). The 



MMRZ is the only model that does not fit the lowest red- 
shifts (2^ 0.3), providing an anti-biased value at the current 
epoch, 6(0) = 0.87. 

• The relative bias difference of the various fitted models 
with respect to that of BPR, A6bpr(z), indicates that the 
BPR, JING and TRK models have a very similar redshift de- 
pendence for z< 2.2 (with |A6 B pr(z)|~ 0.05), while all the 
models show very large such deviations for zlt 2.8, reaching 
|A6bpr|^ 0.3 at the largest redshifts. The SMT and MMRZ 
models show large deviations at the lowest redshifts as well. 

• Beyond the fact that the BPR model provides by far 
the best fit to the QSO bias data, the second best model is 
the TRK model, with AAIC C ~ —2.9. Furthermore, one can 
distinguish that the model pairs (JING, TRK) and (JING, 
MMRZ) are statistically equivalent (AAIC C < 2). 

• The traditional SMT and the recently proposed MMRZ 
models rate the worst among all the other one parameter 
models, but interestingly the former provides consistent val- 
ues of Mh and 6(0) with those of the BPR model. 

We attempt now to provide a robust average value of the 
DM halo mass that hosts optical QSO, using an inverse- AIC C 
weighting of the different one parameter model results. This 
procedure provides a weighted mean and combined weighted 
standard deviation of the DM halo mass of: 

(jiM ht <TM h ) = (2.72,0.56) x 10 12 h^Ms 

while the weighted scatter of the mean is ~ 0.44 x 
10 12 fc _1 Af s . 

Finally, we point out that since it appears that mainly 
the 2 high-z bias points are the ones that give the advantage 
to the BPR model with respect to the others, we perform 
a more conservative comparison among the models by ex- 
cluding these two high-z data points. We find that although 
the resulting halo mass and 6(0) are very similar to those 
of Table 1, with variations of a few percent, there are now 
three models that perform equivalently well, the BPR, JING 
and TRK with AIC C — 13. The other two models perform 
moderately (SMT) or significantly (MMRZ) worse, as was 
the case also in the full data comparison, with AAIC C ~ 2 
and 4, respectively. 

4-1-2 Two free parameter models 

We now allow a second parameter to be fitted simultaneously 
with the DM halo mass. Since the free parameters of the bias 
models have been determined using N-body simulations, it 
would be interesting to investigate if their simulation-based 
value can be reproduced by real observational data. The sec- 
ond free parameter that we will use is ax, A, a, y and a for 
the BPR-I, BPR-II, SMT-I, TRK-I and MMRZ-I models, 
respectively. Note that in the case of the BPR-II model we 
will use the simulation based value of ai, with the free pa- 
rameter A being the halo-merging parameter of the BPR 
model (defined in section 3.5). 

Table 2 presents the best fit model parameters result- 
ing from the x 2 minimization procedure, with the first 
and second columns representing respectively the result- 
ing halo mass, Mh, and the second free parameter, while 
the third column the value of the bias at z = 0. In Fig- 
ure 2 we compare the resulting bias evolution models with 
the WMAP7 scaled QSO bias data (as in Figure 1), while 
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Model 


W 12 h- 1 M e 


2 nd param. 


6(0) 


X 2 ■ /df 


AIC C 


BPR-I 


2.2 ±0.4 


4.64 ±0.07 


1.08 


11.95/20 


16.16 


BPR-II 


2.8 ±0.6 


0.008 ± 0.005 


1.02 


12.45/20 


16.66 


SMT-I 


27.0 ±4.0 


0.40 ±0.02 


0.95 


18.72/20 


22.93 


TRK-I 


0.4 ±0.1 


6.14 ±0.43 


1.11 


13.37/20 


17.58 


MMRZ-I 


0.3 ±0.1 


1.21 ±0.04 


0.87 


21.27/20 


25.48 



Table 2. Results of the x 2 minimization procedure between the 
optical QSO data and the bias models with 2 free parameters. 



in the lower panel we present the relative difference be- 
tween the BPR-II model and all the rest, ie., A6bpr-ii(z) = 
[bi(z) - &bpr-ii(z)]/&bpr-ii(«). 

Below we list the main conclusions of the above fitting 
procedure: 

• A first important result is that the only model that re- 
produces the simulation-based second free parameter value, 
is the BPR-I model. The simulation based value is ai = 4.53 
while the fitted value, based on the QSO bias data, is 
Qi = 4.64 ±0.07. This fact will allow us to derive the depen- 
dence of the parameters of the BPR bias evolution model on 
the relevant cosmological parameters (see Appendix A2). 

• Fitting the BPR-II model to the QSO data provides 
A = 0.008± 0.005 which is almost identical to the simulation 
determined value, used in the BPR case (A — 0.006). As it 
is therefore expected, the fitted values of Mh and 6(0) are 
almost identical to those of the one parameter BPR model, 
but the statistical significance of the BPR-II model is lower 
than that of the BPR due to its 2 free parameters. 

• Although the SMT-I and TRK-I models appear now to 
fit slightly better the QSO bias data, especially the higher z 
range, this happens on the expense of providing unexpected 
values for the Mh and very different values of the second fit- 
ted parameter with respect to their simulation based value. 
For example, the SMT-I model provides a huge halo mass, 
~9 times larger than that provided by the corresponding one 
parameter model. This should be attributed to the fact that 
the fitted second parameter, a, is significantly smaller than 
the nominal value of 0.707. Similarly, the TRK-I model pro- 
vides a very small value of Mh , a factor of ~ 9 less than of 
the corresponding one parameter model, while the resulting 
value y = 6.14 implies that A ~ 10 6 , a value extremely large 
and unphysical. Finally, the MMRZ-I model provides again 
a very small value of Mh, while it is the only two-parameter 
model that fits the data worst than the corresponding one 
parameter model. This is due to the fact that we have used 
k = 0.44 and not k = 0.23, which is used in the one free 
parameter model, as suggested by MMRZ. Had we used the 
latter k value we would have found an extremely small value 
of Mh — 10 10 7i _1 Mq. These results probably indicate a de- 
generacy between the two fitted parameters, a fact which we 
indeed confirm for the SMT-I, TRK-I and MMRZ-I models, 
as can be seen in Figure 3 where we plot the 1, 2 and 3cr 
likelihood contours in the parameter solution plane. Con- 
trary to the above models, no such degeneracy is present 
for BPR-I model. Note that in Fig.3 the cross indicates the 
best 2-parameter solution, while the dashed line indicates 
the simulation based value of the second parameter. In the 
case of the TRK-I model the latter corresponds to the virial- 
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Figure 2. Comparison of the QSO bias data with the two free 
parameter bias models (their line types and colors are indicated in 
the Figure). Lower Panel: The relative difference, A6bpr— n(-z), 
between the BPR-II and the rest of the models. 



ization value (y — log 10 355), used in the one free parameter 
fit. 

• Beyond the previously mentioned fundamental flow of 
the two parameter models (SMT-I, TRK-I and MMRZ- 
I), they all provide relatively comparable to the BPR-II 
model fits of the QSO bias data but only within the range 
0.8;$ 2.4 (see lower panel of Fig. 2). Furthermore, the 
MMRZ-I model, as in the case of the one free parameter 
fit, provides an anti-biased value at 0.2, while it also 
provides the worst overall fit to the QSO bias data. 

• Finally, the BPR one parameter model scores the best 
among all the one or two free parameter models and over the 
whole available QSO bias redshift range, while it is statisti- 
cally equivalent with the BPR-I and BPR-II models (since 
|AAIC C |< 1.8; see Table 3). 

We assess in a more quantitative manner the statistical 
relevance of the different theoretical bias models in repre- 
senting the observational QSO bias data by using the in- 
formation theory parameter AIC C and presenting in Table 
3 the model pair difference AAIC C . As previously discussed 
a smaller AIC C value indicates a model that better fits the 
data, while a small |AAIC C | value (ie., ^ 2) indicates that 
the two comparison models represent the data at a statisti- 
cally equivalent level. Due to the resulting unphysical sec- 
ond free parameter, as discussed previously, we do not use 
in Table 3 a comparison based on the SMT-I, TRK-I and 
MMRZ-I models. It is obvious that the one free parameter 
BPR model fairs the best among any model, while it is sta- 
tistically equivalent, as indicated by the relevant values of 
AAICc, to the BPR-I and BPR-II models and to a slightly 
lesser degree to the TRK model. 
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Figure 3. Contour plots of the two fitted parameter solution 
space. The la level is indicated by the relatively thicker red curve. 
The blue dashed line indicates the simulation based value of the 
y-axis parameter. 



BPR-I 


BPR-II 


SMT 


JING 


TRK 


MMRZ 


BPR -1.3 


-1.8 


-10.3 


-5.1 


-2.9 


-7.3 


BPR-I 


-0.5 


-9.0 


-3.8 


-1.6 


-6.0 


BPR-II 




-8.5 


-3.3 


-1.1 


-5.5 


SMT 






5.2 


7.4 


3.1 


JING 








2.2 


-2.1 


TRK 










-4.3 



Table 3. Results of the pair difference AAIC C for the bias evo- 
lution models fitted to the optical QSO data. 



Model 


10 11 /i- 1 M Q 


6(0) 




AICc 


BPR 


6.0 ± 2.5 


0.99 


0.29/4 


2.29 


SMT 


6.4 ± 1.9 


0.90 


0.45/4 


2.45 


JING 


5.1 ± 1.3 


0.83 


0.93/4 


2.93 


TRK 


7.8 ± 1.9 


0.86 


0.74/4 


2.74 


MMRZ 


6.1 ± 1.3 


0.70 


2.10/4 


4.10 



Table 4. Results of the x 2 minimization procedure between the 
one free parameter models and the optical VVTS galaxy bias 
data. 
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Figure 4. Results of the x 2 minimization procedure between 
the VVTS galaxy bias data and the bias models with one free 
parameter. 



4.2 Optical VVTS galaxy Results 

This model-data comparison takes place at relatively low 
redshifts (z < 1.5), covering a small dynamical range in z, 
and therefore we will use only the one free parameter models 
to fit the galaxy bias data. An additional reason is that even 
with the much larger z-dynamical range covered by the QSO 
data, the second parameter could not be constrained (except 
for the case of the BPR-I and BPR-II models). 

The results of the x 2 - m i mm i za tion procedure are pre- 
sented in Table 4, while in Figure 4 we present the model fits 
to the galaxy bias data. Note that the layout of Figure 4 is 
as Figure 1. It is evident that the BPR model fairs the best 
providing the lowest reduced x 2 and AIC C parameter with 
respect to the other models, while the MMRZ model fairs 
the worst. However, due to the small dynamical range in 
redshift, the information theory pair model characterization 
parameter, AAIC C , indicates that all the bias models are 
statistically equivalent in representing the bias data, since 
AAIC C ^ 1.8 for any model pair. As in the QSO case, we pro- 
vide an average halo mass that hosts VVTS optical galaxies 
using an AIC C weighted procedure over the different one- 
parameter bias models. The resulting weighted mean and 



combined weighted standard deviation are: 

(l*M h ,<TM h ) = (6.3,2.1) x lO^h^M® 



while the weighted scatter of the mean is also ~ 0.9 x 
10 n /i -1 M Q . 

It is interesting to point out that the only model that 
finds that at z = the optical galaxies are unbiased (bo — 1), 
in agreement with other studies of wide-area optical galaxy 
catalogues (Verde et al. 2002; Lahav et al. 2002), is the BPR 
model, while all the other models indicate that optical galax- 
ies are quite anti-biased with 6(0) ^ 0.9. 



5 CONCLUSIONS 

In this work we assess the ability of five recent bias evolu- 
tion models to represent a variety of observational bias data, 
based either on optical QSO or optical galaxies. To this end 
we applied a \ 2 minimization procedure between the ob- 
servational bias data, after rescaling them to the WMAP7 
cosmology, with the model expectations, through which we 
fit the model free parameters. 

In performing this comparison we assume that each halo 
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is populated by one extragalactic mass tracer, being a QSO 
or a galaxy; an assumption which is justified since the obser- 
vational data have been estimated on the basis of either the 
large-scale clustering (> l/i -1 Mpc), corresponding to the 
halo-halo term, or the corrected for non-linear effects vari- 
ance of the smoothed (on 8 Mpc scales) density field. 

The comparison shows that all models fit at an ac- 
ceptable level the QSO data as indicated by the their re- 
duced x 2 values. Using the information theory characteris- 
tic, AlCc, which takes into account the different number of 
model free parameters we find that the model that rates 
the best among all the other is the Basilakos & Plionis 
(200f; 2003) model with the tracer merging extension of 
Basilakos, Plionis & Ragone-Figueroa (2008), which is the 
only model fitting accurately the optical QSO bias data over 
the whole redshift range traced (0 < z < 4) . The only other 
model that is statistically equivalent at an acceptable level 
is that of Tinker et al. (20f0). The average, over the differ- 
ent bias models, DM halo mass that hosts optical QSOs is: 
M h ~ 2.7(±0.6) x lO 12 /i _:L M . 

Finally, all the investigated bias models fit well and 
at a statistically equivalent level the VVTS galaxy bias 
data, with the BPR model scoring again the best, and 
the MMRZ the worst. The average, over the different bias 
models, DM halo mass hosting optical galaxies is: Mh — 
6(±2) x W 11 h~ 1 M @ . 
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APPENDIX A: SIMULATION BASED BPR 
MODEL PARAMETER ESTIMATION 

Al ACDM Simulations 

We have run a new WMAP7 ACDM N-body simulation us- 
ing the GADGET-2 code (Springel 2005) with dark mat- 
ter only. The size of the box simulation is 500/i _1 Mpc 
and the number of particles is 512 3 . The adopted cosmo- 
logical parameters are the following: Q m = 0.273 SI a = 
0.727, h = 0.704, as = 0.81 and the particle mass is 
7.07x 10 10 /i _1 Mq, comparable to the mass of a single galaxy. 
The initial conditions were generated using the GRAFIC2 
package (Bertschinger 2001). We also use a similar size sim- 
ulation, generated in Ragone-Figueroa & Plionis (2007), of 
a ACDM model with Q m = 0.3, S1 A = 0.7, h = 0.72 and 
a s = 0.9. 

The dark matter haloes were defined using a FoF algo- 
rithm with a linking length I = 0.17(?i)~ 1//3 , where (n) is 
the mean particle density. 
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Figure Al. The Parameters Ci and C2 derived from the WMAP7 
(filled points) and WMAP1 (open points) ACDM simulation 
(points). Continuous lines correspond to the function form of eq. 
(A3), with best fit parameters shown in Table Al. 



Model 


on 


a 2 


fa 


fa 


WMAP1 


3.30±0.13 


-0.36±0.01 


0.34±0.04 


0.32±0.04 


WMAP7 


4.53±0.22 


-0.41±0.02 


0.37±0.04 


0.36±0.04 



Table Al. Results of the \ 2 minimization used to evaluate the 
parameters that enter in the Ci and C2 constants (which depend 
on DM halo mass) of the BPR bias evolution model. 
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Figure A2. Correlation between (a) the best fitted c?i values 
of the BPR model using the QSO bias data scaled to different 
cosmologies for a grid of f2 m and erg values and (b) the predicted 
ai values, based on eq, l|A4JI . The red points correspond to Q m ^ 
0.273 while the black points to Q m > 0.273. 



We estimate the bias redshift evolution of the different 
DM haloes, with respect to the underlying matter distribu- 
tion, by measuring their relative fluctuations in spheres of 
radius 8 h~ Mpc, according to the definition of eq.([3J), ie., 



b(M, z) = 



o&,h{M,z) 
»(*) 



(Al) 



where the subscripts h and m denote haloes and the under- 
lying mass, respectively. The values of ag : h{M, z), for haloes 
of mass M, are computed at different redshifts, z, by: 



al h {M,z) = 




(A2) 



where N is the mean number of such haloes in spheres of 8 
h^ 1 Mpc radius and the factor 1/N is the expected Poisso- 
nian contribution to the value of a\ h . Similarly, we estimate 
at each redshift the value of the underlying mass <Ts,m- In 
order to numerically estimate a 2 j we randomly place -/V ran d 
sphere centers in the simulation volume, such that the sum 
of their volumes is equal to ~ 1/8 the simulation volume 
(A'rand = 8000). This is to ensure that we are not oversam- 
pling the available volume, in which case we would have 
been multiply sampling the same halo or mass fluctuations. 
The relevant uncertainties are estimated as the dispersion 
of erf j over 20 bootstrap re-samplings of the corresponding 
halo sample. Note that we do not explicitely correct for 
possible non-linear effects in 8 (although the density field is 
indeed smoothed on linear scales - 8 /i -1 Mpc); we do how- 
ever expect that such effects should be mostly cancelled in 
the overdensity ratio definition of the bias. 

We use the DM halo bias evolution, measured in the two 
simulations, for different DM halo mass range subsamples in 
order to constrain the constants of our bias evolution model, 
ie., Ci, C2. The procedure used is based on a \ 2 minimization 



of whose details are presented in BPR and thus will not be 
repeated here. In Fig.Al we present as points the simulations 
based values of these parameters, for both cosmologies used, 
and as continuous curves their analytic fits, which are given 
in eqs. (|20p and (|2ip . The resulting values of the parameters 
01,2 and /3i,2 can be found in Table Al. It is interesting to 
note that the slope of the functions Ci and Ci is roughly a 
constant and independent of cosmology, with a value fa ~ 
fa ~ 0.35(±0.06). 

A2 Dependence of the BPR model constants on 
Cosmology 

The dependence of the constants ot\ and 02 of the BPR 
bias model on the different cosmological parameters is an 
important prerequisite for the versatile use of the model in 
investigating the bias evolution of different mass tracers and 
determine the mass of the dark matter halos which they in- 
habit. In Basilakos & Plionis (2001) we predicted a power 
law dependence of Q2 on fl m . Indeed, fitting such a depen- 
dence, using the WMAP1 and WMAP7 ACDM simulations, 
we find: 

/0 273 \ " 

a 2 (n m ) = -0.41 (— — with n~ 2.8/2 (A3) 

consistent with the value n — 3/2 anticipated in Basilakos 
& Plionis (2001). 

Now, in order to investigate the dependence on different 
cosmological parameters of the parameter ai, we have used 
the optical QSO data and the procedure outlined in section 
2.2 to scale the QSO bias data to different flat cosmologies, 
using a grid of fl m and erg values. The grid was defined as 
follows: n m G [0.18,0.5] and cr 8 G [0.7,0.94], both in steps 
of 0.01. We then minimize the BPR bias evolution model 
to the scaled bias data to different cosmologies QSO, finally 
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providing for each pair of (Q m , fs) values the best fitted «i 
and Mh values. 

Then using a trial and error approach to select the 
best functional dependence of the derived Qi(f2 m ,crg) val- 
ues to the relevant cosmological parameters, we find a best 
fit model of the form: 



In Fig. A2 we correlate the derived ai(f2 m , as) values, result- 
ing from fitting the BPR bias evolution model to the scaled 
QSO bias data, to those predicted by the model of eg . (I A ll) . 
It is evident that the correspondence is excellent, indicating 
that indeed the above estimated cosmological dependence of 
ai is the indicated one. 




(A4) 



with 




Q m 0.273 
n m > 0.273 



(A5) 
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